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Many real phenomena may be modelled as random closed sets in W^, of different Hausdorff di- 
mensions. In many real applications, such as fiber processes and 7i-facets of random tessellations 
of dimension n < d in spaces of dimension d > 1, several problems are related to the estimation 
of such mean densities. In order to confront such problems in the general setting of spatially 
inhomogeneous processes, we suggest and analyze an approximation of mean densities for suf- 
ficiently regular random closed sets. We show how some known results in literature follow as 
particular cases. A series of examples throughout the paper are provided to illustrate various 
relevant situations. 
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1. Introduction 

Many real phenomena may be modelled as random closed sets in R'', of different Hausdorff 
dimensions (see, e.g., [3, 4, 13, 17, 19, 23, 24]). 

We recall that a random closed set Q in is a measurable map 

e:(f},J-,P)^(F,aF), 

where F denotes the class of closed subsets in and crjr is the tr-algcbra generated by 
the so-called hit-or-miss topology (see [16]). 

Let 0„ be almost surely (a.s.) a set of locally finite Hausdorff n-dimensional measure 
and denote by the n-dimensional Hausdorff measure on M.'^. The set 0„ induces a 
random measure /ie„ defined by 

Me„(A):=H"(e„nA), Ae 
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(for a discussion of the delicate issue of measurability of the random variables ?^"(0„ n A), 
we refer to [2, 26]). 

Under suitable regularity assumptions on a random closed set 6„ C R'^ with locally 
finite n-dimensional measure a.s., in [7] the concept of mean density of 0„, that is, the 
density of the expected measure 

E[MeJ(^):=E[7^"(e„nA)], AeB^., (1) 

with respect to the standard Lebesgue measure i/'^ on M.'^, has been revisited in terms of 
expected values of a suitable class of linear functionals (Delta functions a la Dirac) . 

It is clear that if n < d and /xe„(w) is a Radon measure for almost every uj Gfl, then it 
is singular with respect to the d-dimensional Lebesgue measure v"^. On the other hand, in 
dependence of the probability law of 0„, the expected measure may be either singular or 
absolutely continuous with respect to ly'^; in the latter case, we say that a random closed 
set is absolutely continuous in mean, and its mean density is the classical Radon-Nikodym 
derivative of E[At0„] with respect to v'^. 

In many real applications, such as fiber processes and n-faccts of random tessellations 
of dimension n < d in spaces of dimension d > 1, several problems are related to the 
estimation of their mean densities (see, e.g., [3, 5, 12, 18, 19, 24]). The scope of this 
paper is to provide concrete approximations that may lead to (possibly local) estimators 
of mean densities of lower-dimensional inhomogeneous random closed sets; the globally 
stationary case has been the subject of various papers and books in the literature (see, 
e.g., [3]), while the problem of local density estimation for n = (absolutely continuous 
random variables) has been largely solved for quite some time - this now being part of the 
standard literature - by means of either histograms, or kernel estimators (see, e.g., [10, 
22]). In order to confront such problems in the general setting of spatially inhomogeneous 
processes, we suggest and analyze here an approximation of mean densities for sufficiently 
regular random closed sets Qn in R''- We shall show how some known results in the 
literature follow as particular cases. 

A series of examples throughout the paper are provided to illustrate various relevant 
situations. 

2. Preliminaries and notation 

In this section, we collect some basic facts and terminology that will be useful in the 
sequel. 

We shall define a Radon measure in to be any non-negative and cr-additive set 
function /i defined on the Borel cr-algebra B-ig^d which is finite on bounded sets. 
We know that every Radon measure /i on M.'^ can be represented in the form 

where fj,<^ and fj,± are the absolutely continuous part of fi with respect to i/'^ and the 
singular part of fj,, respectively. Denoting by Br{x) the closed ball with center x and 
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radius r, as a consequence of the Besicovitch derivation theorem (see [1], page 54), we 
have that the hmit 

exists in M for v'^-a..c. x £ R'' and it is a version of the Radon-Nikodym derivative of /^<g, 
while /i^ is the restriction of ^ to the z/''-neghgible set {a; G R*^ : Sfi{x) = cxd}. We caU 
the density of the measure /i with respect to ly"^. 

Let us consider a random closed set in M.'^ with locally finite n-dimcnsional measure 
a.s., say 0„, such that the expected measure E[/^e„] induced by 0„ defined by (1) is a 
Radon measure and let (5E[/ie„] be its density. We refer to (5E[^e„] as the mean density 
of 6n- (For a more detailed discussion about the density of deterministic and random 
closed sets, see [7].) For any random closed set 0„ in R"* with locally finite n-dimensional 
measure for some n < d a.s., while it is clear that /^e„(w) is a singular measure, it may well 
happen (e.g., when the process Qn is stationary) that the expected measure E[/i0^] is 
absolutely continuous with respect to i^'', so that (5E[//e„] is the Radon-Nikodym density 
of E[/ie„]- In this case, the following definition is given [6]. 

Definition 1 (Absolute continuity in mean). Let 8„ he a random closed set in 
with locally finite n-dimensional measure a.s. such that its associated expected measure 
E[/ie„] is a Radon measure. We say that 0„ is absolutely continuous in mean i/E[/ie„] 
is absolutely continuous with respect to v'^ . In this case, we denote by Ae„ its density. 

Remark 2 (The Q- dimensional and d-dimensional cases). If n = and 8o(w) = 
{X{uj)} is a random point, then E[H"{eo n •)] = G •)• Therefore, So is absolutely 
continuous in mean if and only if the law of X is absolutely continuous and Aeo coincides 
with the probability density function of X . 

On the other hand, E[/ied] is always absolutely continuous with respect to t^'* and by 
a simple application of Fubini's theorem (in f2 x R^, with the product measure P x ly'^), 
we have that [21] (see also [15]) 



E[/xe, (B)] = / P{x G Grf) d.T VB G 

JB 



SO Ae^ {x) = P(.T G Qd) for v'^-a.e. cc G R''. 

We mention that a stronger definition of absolute continuity in mean of a random 
closed set has been given in [6] in terms of the expected measure of its boundary. 

The aim of the present paper is to provide an approximation of the expected measure 
E[/.Je„] under quite general regularity assumptions on the random closed set 0„, which 
provides, as a by-product, an approximation of the mean density Ae„ when E[^e„] is 
absolutely continuous with respect to v^^ . To this end, the notion of Minkowski content 
of closed sets in R** plays here a fundamental role. Denoting by bk the volume of the unit 
ball in R'^ and by S^r the closed r- neighborhood of S*, that is. 

S'er := {x G M'' : 3?/ G S" with \x~y\< r} 
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(known also as the Minkowski enlargement of S with the closed ball i?,.(0)), the n- 
dimensional Minkowski content of a closed set S C Mf^ is defined by 



r4,0 bd-nr'^-" ' 

whenever the limit exists. 

General results concerning the existence of the Minkowski content of closed subsets in 
M'^ are well known in the literature, related to rectifiability properties of the sets involved. 
Given an integer < n < d, we say that a set C C K'^ is countably T-V^ -rectifiable if there 
exist countably many Lipschitz maps fi : M" — ^ M'' such that 

H"|^c\y =0. 

Rectifiable sets include piecewise sets, and they still have nice properties from the 
measure-theoretic viewpoint (e.g., one can define an n-dimensional tangent space to them, 
in an approximate sense); we refer to [1] for the basic properties of this class of sets. We 
quote the following result from [1], page 110. 

Theorem 3. Let S CR.'^ be a countably -rectifiable compact set and assume that 

7]{Br{x))>jr'' VxG5,VrG(0,l) (2) 

holds for some 7 > and some Radon measure rj in absolutely continuous with respect 
to H". Then 

limZ!%L^H"(5). 
In particular, in Lemma 6, we shall prove a local version of the above theorem. 



3. Approximation of mean densities 



Throughout the paper, we consider countably ?^"-rcctifiable random closed sets in 
(i.e., P-a.e. w G il, 6n(i^) C M'' is a countably 7{"-rectifiable closed set), with Radon 
measure E[/ie„]. 

By the aforementioned Besicovitch derivation theorem, a natural approximation at the 
scale r of Ae„ can be given by 



bdr^ 



It can be easily proven that the associated measures J^'''E[/^e„]i''^ weakly-* converge to 
the measure E[//e„] as r 4, (i.e., in the duality with continuous and compactly supported 
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functions). This convergence result can also be understood by noting that ^^''■'E[/xe„](x) 
is the convolution of the measure E[/^e„] with the kernels (here, 1e stands for the char- 
acteristic function of E) 

Priy) ■■= j^^BAo)iy)- 

Noting that a point y belongs to Br{x) if and only if the point x belongs to the 
Minkowski enlargement y^r of it follows that whenever is a random point in , 

E[H°(eoni3.(x))] meo^M)] 
5^ ^E[Meo](x) = b.r'^ ■ 

In the particular case d = 1 with Go = X a, random variable, we have that 

5W]E[^^](,).ffi£tli:^; (3) 

if X is a random variable with absolutely continuous law and pdf fx, we have (see Remark 
2) that Ax = fx, so that (3) leads to the usual histogram estimation of probability 
densities (see [20], Section VII. 13). 

We point out that even if a natural approximation of E[/xe„] can be given by the func- 
tions (5('')E[/ie„] defined above, problems might arise in the estimation of E[^e„ (^r(a;))], 
as the computation of the Hausdorff measure ■H"(0„(a;) r\Br{x)) is typically non-trivial. 
As a matter of fact, a computer graphics representation of lower-dimensional sets in 
is anyway provided in terms of pixels, which can only offer a 2-D box approximation 
of points in (an interesting discussion on this is contained in [14]). Therefore, we 
are led to consider a new approximation, based on the Lebesgue measure (much more 
robust and computable) of the Minkowski enlargement of the random set in question. 
This procedure is obviously consistent with (3). We now state the main theorem of this 
paper. 

Given a random closed set 0„ with locally finite n-dimensional measure a.s. in M'* (the 
particular case n = d is trivial), for any compact window W C M'', let Twi^n) '■ ^ — > ^ 
be the function defined as follows: 

rw(0ri) ■ = sup{7 > : 3 a probability measure rj <^ such that 

(4) 

v{Br{x)) > 7r" Vx- G e„ n VFei,Vr e (0, 1)}. 

Theorem 4 (Main result). Let 6„ be a countably H" -rectifiable random closed set in 
M'' such that E[/ie„] is a Radon measure. Assume that for any compact window W C R'', 
there exists a random variable Y with E[F] < oo, such that l/Vw(Qn) < Y a.s. Then 

for any bounded Borel set A C M'' such that 

E[^eJ(aA) = 0. (6) 
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In particular, if 0„ is absolutely continuous in mean, we have 

lim / ^(^gQ;e^) f {x)dx (7) 

for any bounded Borel set A C with v'^{dA) = 0, where Ae„ is the mean density of 
Qn ■ Finally, if 0„ is stationary, we have 

j.^P(xoee^^ VxoeM^ (8) 

having denoted by Lq^ the constant value of Xq^ . 

Remark 5. In the statement of the above theorem, we introduced the auxihary random 
variable Y in order to avoid the non-trivial issue of the measurability of Tw; as a matter 
of fact, in all examples, one can estimate l/Fw from above in a measurable way. 

Note also that if 0„ satisfies the assumption of the theorem for some closed W , then 
it satisfies the assumption for all closed W C W; analogously, any random closed set Q'^ 
contained almost surely in 0„ still satisfies the assumption of the theorem. 

Finally, note that condition (6), when restricted to bounded open sets A, is "generically 
satisfied" in the following sense: given any family of bounded open sets {At}teA with 
closes C At for s < t (clos A denotes here the closure of the set A), the set 

T:=={tGR:E[/xeJ(5At) >0} 

is at most countable. This is due to the fact that the sets {dAt}t^T are pairwise 
disjoint and all with strictly positive E[/ie„]-measure (in particular, the sets {t G 
(— cxd,to) :E[/ie„](9At) > l/m} have cardinality at most mE[/ie„](^m))- 

3.1. Proof of the main result 

Note that by a straightforward application of Fubini's theorem, the limit in (5) is equiv- 
alent to 

which might be regarded as the local mean n-dimensional Minkowski content of 0„. 
The problem hence reduces to finding conditions on 6„ which ensure that (9) holds. 
To this end, let us observe that if Qn is such that almost every realization 0„(aj) has 
n-dimensional Minkowski content equal to the n-dimensional Hausdorff measure, that is. 



(10) 
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then it is clear that, taking the expected values on both sides, (9) is strictly related to the 
possibility of exchanging limit and expectation. We therefore ask whether (10) implies a 
similar result when we consider the intersection of 0„^^(cj) with an open set A in R'', 
and for which kind of random closed sets the convergence above is dominated, so that 
exchanging limit and expectation is allowed. 

The following result is a local version of Theorem 3. 

Lemma 6. Let S be a compact subset of M'* satisfying the hypotheses of Theorem 3. 
Then, for any A 6 B-g^d such that 



the following holds: 



7^"(5naA) = o, (11) 



^.^ ^^(5e nA) ^^^ (12) 



Proof. If n~d, then equality (12) is easily verified. Thus, let n < d. 

We may note that by the definition of rectifiability, if C C M'* is closed, then the 
compact set S CiC is still countably H"-rectifiable; besides, (2) holds for all points x S 
S nC (since it holds for any point x G S). As a consequence, by Theorem 3, we may 
claim that for any closed subset C of R'', the following holds: 

i^ ^'iisr^c)^r) (13) 

no bd-nr'^-"^ 
Let A be as in the assumption. 

• Let £ > be fixed. We may observe that the following holds: 

S'er n A C (S'nclosA)0r U (S'nclosAee \int A)0r Vr < e, 

where int A denotes the interior of the set A. 

Indeed, if a; £ S,^r H A, then there exists y E S with \x — y \ < r and y e clos ^0^. If 
X ^ {S O closA)0r, we must then have y E S \ closA, hence y € S D clos \ closA. 

By (13), since closA and closA^^ \intA arc closed, we have 

lim -'(iS^closA)^.) ^ ^ ^^^^^^ (IJ) ^ ^^^^ 

I^'^(S'nclosA0£ \ intA)0r ^,ri , , A X ■ A^ 

lim — ^ 7^"(S'nclosA0e \int A). (15) 

r\.0 Od—nT 



Thus, 



v''{S^rr\A) 

hm sup — -z 
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v'^US nc\osA)^r U (^ndosAee \ int A)er) 
< lini sup 



z^'^((S'nclosA)©r) + i^''((<S'nclosAee \iiitA)er) 

p 

r4,0 

(14).(15). 



< lim sup ■ , . 



•H" (5 n A) + -H" (S" n clos Aee \ int A) ; 
by taking the limit as e tends to 0, we obtain 

v'^is^^nA) 



lim sup ^ < -H" (5 n A) + 7^" (S* n 9/1) = (5 n A) . 



=0 

• Now, let i? be a closed set well contained in A, that is, the Hausdorff distance 
between A and B is greater than 0. There then exists f > such that S^,, C A, V?' < f. 
So, 

H"(5nS)^L3himinf4(^I2#M 

< lim mt — ; 3 

< liminf^^^(^-fl^. 

Let us consider an increasing sequence of closed sets {-Bm}mGN well contained in A such 
that Bm /^intA. By taking the limit as m tends to oo, we obtain that 

liminf ^^(^g'^l^^) > lijn (5 n B^) = n int A) ^ = ' (5 n A). 



We summarize: 

^ '-msu" 

and so the thesis follows. '-' 



(5 n A) < lim inf < lim sup ^y^^^''^' < (5 n A) 



If we consider the family of random variables h^^^-n. ^ ; for r going to 0, we ask 
which conditions must be satisfied by a random set 0„ so that they are dominated by 
an integrablc random variable. In this way, we could apply the dominated convergence 
theorem in order to exchange limit and expectation in (12). 

Lemma 7. Let K be a compact subset o/R'^ and assume that 

7?(B^(x)) >7r" Vxe A',Vr e (0,1) (16) 
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holds for some 7 > and some probability measure 77 in R**. 
Then, for all r <2, 

Proof. Since /^©r is compact, it is possible to cover it with a finite number p of closed 
balls -83^(0;^), with Xi G -ftr©r, such that 

\xi-Xj\>3r, i^j. 

As a consequence, there exist yi,. . . ,yp such that: 

• yi€K,i^l,...,p; 

• \yi-yj\ >r, i^i; 

• KfSr^{f^=lBir{y^)■ 

In fact, if Xi G if , then we choose yi ~ Xi ; if Xi e K^r \ K , then we choose yi G Br {xi) C\K . 
As a consequence, \yi — Xi\ <r and B^riyi) 12 B^ri^i) for any i ~ I, . . . ,p. So, 

p p 

[jBiriy^)^[jBsriXi)^Ks,r 
i=l i=l 

and 

3r < \xi -Xj\< \xi -yi\ + \yi -yj\ + \yj -X]\< 2r + -yj\, i^j. 

For r < 2, -6^/2 (j/i) H Bj./2iyj) = <ZS. Since, by hypothesis, 77 is a probability measure 
satisfying (16), we have that 



/ P \ ^ (16) 

l>v([JBr/2{y^)] =J2'^iBr/2iy^)) > PI 
\i=l I i=l 

and so 

1 2" 

ry ipil' 



In conclusion, 



bd-nr'^-'' - bd-nr'^-'" " " 7 bd-n ' ° 

We are now ready to prove the main result of the paper. 

Proof of Theorem 4. Let be a compact window in M.'^ and Twi'dn) the function 
defined in (4). Let Y be as in the assumptions of the theorem and A be any Borel subset 
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of R'^ such that 

AcintW and E[-H"(e„ n = 0. (17) 

Smce E[Y] < oo, we have Y{ijj) < oo for P-a.e. u eil. With no loss of generahty, we shaU 
assume that l/Twi^n) < Y a.s. 
Let us define 

Ha := G f7 : -H"(e„(cj) n dA) = 0}, 

VIt := {cl> e : 0„(cl') is countably ?^"-rectifiable and closed}, 
fiy {cj e : Y{ijj) < oo}, 

f^r:= l^en:—-^-—<Y{uj) 

by hypothesis, P(17a) = P(Ot) P(r2y) = P(r2r) = 1- 

Thus, if n' := 17^ n n rjy n fir, it follows that = 1. 

Let cl> e ri' be fixed. Then: 

• 1/Y{uj) < rvF(0n(w)), that is, a probability measure 77 ^ T^'' exists such that 

viBrix)) > T7^(c^)r" Vx e e„(^) n W^i^r e (O, l); 
r [uj) 

• H"(e„(cj)n(9^) = 0. 

So, by applying Lemma 6 to 9„ n W(^i and noting that (0„ n M^ei)©r n A = 0,1©^ n A, 
we get 

>Aen^)nA) 

r'iO bd-nr'^ 

that is, we may claim that 



% \ d-n = n A); 



lim4^^) 



lim ' — ^ = 7^" (e„ n A) almost surely. 



Further, for all lj G O', 6ri(w) n W^i satisfies the hypotheses of Lemma 7 and so 

iy^{e,-,^^{uj) n A) i.'^((e„(c^) n i^eOer n ^) 



Let Z be the random variable defined as follows: 

Z{io):^Y{io)2"4'^-^, a; G 17'. 
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By assumption, E[Z] < oo, so that the dominated convergence theorem gives Equation 
(9), and hence the first statement of the theorem. 

The second statement is a direct consequence of the first, because if 8„ is absolutely 
continuous in mean, by Definition 1, we have that E[yUe„](5A) = if u'^{dA) = and 

E[MeJ(A)- / XeAx)dx. 

J A 

Finally, if 0„ is stationary, then P(x € ©n^^) is independent of x and the expected 
measure E[/ue„] is translation invariant, that is, Ae„(a;) = ie„ S K+ for i/'^-a.e. x G K*^. 
It follows that 

lim / ^^-^>)d.^lim^(-°^>).^(A) 
for any xq G R'' and so (8) follows directly by (7). □ 

A classical criterion (see, e.g., [11] or [1]) states that fin weakly-* converge to ^ if and 
only if fin{A) — > fJ-{A) for any bounded open set A with fJ-{dA) = 0. Hence, for any random 
closed set 0„ with locally finite rt-dimensional measure a.s. satisfying the assumptions 
of Theorem 4, it follows that the measures fi®^ defined by 

JA ^d—n' 

weakly-* converge to the measure E[/ie„] as r J, 0. As a consequence, the associated 
density functions 

can be taken to constitute a weak approximating family of the density (5E[/xe„] of 0„ 
when the process is absolutely continuous in mean. 

Remark 8 (Mean density as a pointwise limit). It is tempting to try to exchange 
limit and integral in (7), to obtain 

lim— -r^ — = Ae„(a;), (18) 

at least for i/'^-a.e. x € 'R'^. The proof of the validity of this formula for absolutely contin- 
uous (in mean) processes seems to be quite a delicate problem, with the only exception 
being stationary processes; this is the main reason why the existing literature has exten- 
sively considered this case [3]; we nonetheless wish to stress that global stationarity is 
only a sufficient condition for (18) to hold. Indeed, we already know that in the extreme 
cases n = d and n = 0, it is not hard to prove it. 

In the case n = d, we know from Remark 2 that \q^{x) = P{x G Qd) for ly'^-'A.e. x, and 
obviously F{x G ©d©?.) converges to P(a; G 9^) for all x. 
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In the case n = 0, let 80 = {X}, with X an absolutely continuous random point in M'* 
with pdf fx , and note that 



(xeOoer) nxeBr{x))_ 1 



bdr'^ bd.r'^ bdr 



,d 



fxiy)dy. 

Br(x) 



Therefore, (18) with Xq^ = fx holds at any Lebcsgue point of fx (i-C-, a- point x such that 
the mean values of ?/ 1— > \fx{y) — fx{x)\ on Br{x) are infinitesimal as r ],0) and hence 
for v'^-a..e. xGWK This is the main reason why the problem of local density estimation 
has already been solved [10, 22]. 

We may note that 

P(:r e On^,.) = P(e„ n Br{x) ^0)= Te,. 

thus making explicit the connection to Te„ , the capacity functional characterizing the 
probability law of the random set 0„ [16]. Therefore, the family of functions (5®'' may 
suggest estimators of the mean density of 0„ in terms of the empirical capacity functional 

of e„. 

For instance, even if the special case n = can be handled with much more elementary 
tools, we may note how the well-known and well-studied case of a random point X in M.'^ 
is consistent with our framework. It is immediate to check that X satisfies the hypotheses 
of Theorem 4 (it is sufficient to choose 77 7i'^{X{uj) D ■)) and that, in the particular 
case d~l with X absolutely continuous r.v. with pdf fx, (18) becomes 

P{X e[x-r,x + r]) ^ 
lim 7^ = fx , 

which leads to the usual probability density estimation by histograms. 



4. Applications 

In many real applications, 0„ is given by a random collection of geometrical objects, so 
that it may be described as the union of a family of random closed sets Ei in R'^, 

On^ljE,, (19) 

i 

and a problem of interest is to determine or estimate the mean density Ae„ of Q„ - Such 
unions of random sets can be taken as models for the so-called particle processes, widely 
studied in stochastic geometry. 

We may note that Theorem 4 seems to require sufficient regularity of the i?i's, rather 
than specific assumptions about their probability law or the stochastic dependence among 
them. In the present section, we provide concrete examples and results in this direction, 
in order to clarify the wide range of applicability of our main theorem and to show how 
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some known results in the current literature for stationary random closed sets follow 
from (8). 

For the sake of simplicity, we start by considering the case in which the union (19) of 
the Ei is finite a.s. so that there exists a positive integer-valued random variable $ such 
that 

4> 

e„ = [jE,. (20) 

4=1 

In the next example, we consider a class of random sets of this kind, which might be taken 
to model a class of birth- and- growth stochastic processes, relevant in real applications, as 
discussed in Example 2. 



Example 1 . A class of random sets satisfying the hypotheses of Theorem 4 is given by 
all sets Qn which are a random union of countably 7^"-rectifiablc random closed sets in 
W^, as in (20), such that 

(i) E[$] < oo; 

(ii) Ei,E2,. ■ ■ are i.i.d. as E and independent of $; 

(iii) E[H"{E)] = C < oo and 37 > such that for any wen, 



7^"(£;(w)nB,.(x)) >7r" V.t e £:(w), Vr G (0, 1). 



(21) 



We can choose 77(-) := ^I^q^"^)')-' fixed w e il. As a consequence, 77 is a prob- 

ability measure absolutely continuous with respect to "H" and such that 



rjiBr{x))> 



7 



■H"(e(w)) 



yxeOniuj), re (0,1). 



In fact, if X G 0„(a;), then there exists an 1 such that x G E!i{lu); since 0„(aj) 
{jfli^ £'i(w), we have 

W'{en{w)nBr{x)) ^ W^{E,{uj)nBr{x)) ^ 

r]{Br{x)) = ——: -^^ > . .. > 



H"(e„(c^)) 



H"(e„(w)) 



•H"(e(w)) 



As a resuh, the function r(e„) defined in (4) is such that r(e„) > j/K"{e„) 1/Y 
and so it remains to verify only that E[H"(8„)] < 00: 



E[?^"(e„)] - E[E[7^"(e„)|$]] 



fe=l 



n"i[jE, 



\i=l 



P($ = fc) 



(ii) 

< ^^E[H"(£;,)]P($ = fc) 

fc=l i=l 
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fc=i 

= CE[$] < oo. 

Note that wc have not made any particular assumption on the probabihty laws of $ 
and E. Further, it is clear that the same proof holds, even in the case in which the Ei's 
are not i.i.d., provided that E[H"(-Ei)] < C for all i and (21) is true for any Ei (with 7 
independent of w and i). 

In the next example, we show how unions of time-dependent random closed sets may 
be used to model a class of time-dependent geometric processes as well; it also provides 
an example of an absolutely continuous in mean, but not stationary, random closed set. 

Example 2 (Birth- and- growth process). Roughly speaking, a hirth- and- growth pro- 
cess is a family {6*}t of random closed sets which develop in time, according to a given 
growth model, from points {nuclei^ or germs) that are born at random both in space and 
time (for details, see [6] and references therein). These kinds of processes are described 
by dynamic germ-grain models (see, e.g., [9, 24]), that is, once born, each germ gener- 
ates a grain subject to surface growth with a speed G which is, in general, space-time 
dependent. The nucleation process {T„, X„}„gN, where r„ is the R_)_-valued random vari- 
able representing the time of birth of the nth nucleus and X„ is the M'^-valued random 
variable representing the spatial location of the nucleus born at time Tn, is commonly 
modelled by a marked point process A^. In particular, iV is a random counting measure 
on R-|_ X such that, for any fixed t € R+, the random number Af([0,t] x R'') of nuclei 
born up to time t is finite with probability 1. 

For the sake of simplicity, we assume here that N is given by an inhomogeneous Poisson 
point process in R_|_ x R"^ with intensity a(i,x) (i.e., E[iV(d(t, x)] = a{t,x) dtdx) and that 
the growth occurs with a constant normal velocity G > so that, for any fixed time t, 
8* is a finite union of random balls in W^: 

©*= U ^G(t-T.)(^.)- 
i:Ti<t 

Considering now the random closed set 90*, we observe that it is absolutely continuous 
in mean, but not stationary, for any t > Q. We may prove this (by contradiction) as 
follows. 

Assume E[?^'^~^(99* n •)] to be not absolutely continuous with respect to v^^; there 
then exists some AcW^ with v'^{A) = such that E[H''~H50* ^ A)] > 0. 
It is clear that 

EiW^-^dQ* n A)] > ^ PiW^-^dQ' n A) > 0) > 

and 

r{n'^-\de* n A) > 0) <p{3{Tj,Xj) ■.■H'^-\dBGit-T,){Xj) n A) > 0). 
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As a consequence, we have 

E[n''-\de* n A)] > p($(^) 7^ o) > o, 

where 

A := {{s, y) e [0, <] X M-^ : H^"\dBGit-s) iv) n A) > 0}. 

Denoting by As ■= {y £ M** : (s, y) e A} the section of A at time s and by A^ :— {s > 0: 
(s,y) G y^} the section of A at y, we note that ^^(^'') = for all y because v'^{A) = 
(it suffices to use spherical coordinates centered at y to obtain that the v^-a..e. ball with 
radius s centered at y intersects A in an H''" ^-negligible set). Therefore, we may apply 
Fubini's theorem to get 

/•oo /"OO f p poo P 

/ v\As)As= / XA.dyds= / XAydsdy= (^^) dj/ = 0. 

Jo Jo JTBL'' JM'' Jo Jm'' 

It follows that i'^{As) = for i^^-a.s. s G [0,t] and so 

E[3>(^)] = [ a{s,y)dsdy^ [ [ a{s,y)dyds^O. 
J A Jo J A, 

But this is an absurd, since 

P($(^) 7^ 0) > E[^{A)] > 0. 

Observing that Theorem 4 applies to 90* (see [25], Chapter 4 for details), as a conse- 
quence of (7), we have 

lim / — ZMlI ^ [ Age* i^) dx 

rio 2r J A 

for any bounded Borel set A C M'^ with v''-{dA) = 0, where A^e* is also known as the 
mean surface density associated with the birth-and-growth process Mean volume 

and surface densities of birth-and-growth processes are of a great interest in applications 
(see [6] and references therein). 



Coming back to consider the basic model described in Example 1, we now prove a 
result which relates the probability that a point x belongs to the set Qn^^- to the mean 
number of Ei which intersect the ball Br{x), and which might therefore be useful in 
statistical applications. 

Proposition 9. Let n < d, let ^ be a positive integer-valued random variable with E[$] < 
oo and let {Ei} be a collection of random closed sets with locally finite n-dimensional 
measure a.s. Let 9„ be the random closed set defined as 

e„ = U ^- 

i=l 
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If El, E2, ■ ■ ■ are i.i.d. as E and independent of $, then, for any x G such that P(a: £ 
E)^Q, 

li- ^-P^ = ^^^f = m] (22) 
provided that at least one of the above limits exists. 



Proof. The following chain of equalities holds: 

/ * \ 



^l-¥{[]{xiE,^^}^ 

00 / k \ 

fe=i \i=i / 
00 

= l-^[P(x^£;e.)]'P($ = fc), 



fe=i 



the ii'i's being i.i.d. and independent of $. Denoting by G the probability generating 
function of the random variable $ and s{r) := P(.t G -E'er), we get that 



Xx-ee„^J = l-G(l-s(r)). 



(23) 



We now observe that 



n#{E^■.x&E,^^}] = Y,^ 



k=l 



P($ = fc) 



= ^fcP(.TeE;er)P($ = fc) 



(24) 



fc=i 



= E[$]P(a; e E;er) 



and recall that E[$] = G'(l) and 1 = G(l). By hypothesis, we know that s{r) J, as r J, 0. 
Thus, by (23) and (24), we have that if s(r) = for some r > 0, then P(a; S Onen) = 
¥J[^{Ei:x G Ei^j^}] = for all < R and so the assertion trivially follows, whereas if 
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s{r) > for all r > 0, we have 



lim 



~ lim 



G(l)-G(l-s(r)) 



rioE[#{Er.x€E,^^}] rio G'(l)s(r) 



1 .ii^G(l-sW)-G(l) 



-sir) 



In conclusion, we obtain 
lim 



lim 



G'(l) riO 

ni^{E,:xeE,^^}] 



(25) 



■40 E[#{E,:xeE,^J] 
(25) ^^^^EmE,:xeE,^^}] 



d—n 



^ = ' E[$] lim 



(26) 

(27) 
□ 



Note that for any random closed set 0„ as in Proposition (9), by (22), we infer that 
the probability that a point x belongs to the intersection of two or more enlarged sets 
Ei is infinitesimally faster than r'^"". In fact, denoting by Fr{x) this event and by 
XA ■ ^ {0, 1} the characteristic function of an event A, we have 



so that, since (24) gives 



E 



= E[$]P(.TG-B®.), 



by taking expectations on both sides and dividing by r'' wc get P{Fr{x))/r'^ " — > 
as r 10. 

Remark 10. Whenever 0„ is absolutely continuous in mean and it is possible to ex- 
change limit and integral in (7), we can use the fact that A is arbitrary to obtain 

Ae„ {x) = E[$] hm ^^^^ = E[$] {x) 

for z^'^-a.e. x G R'', where Ae„ and A^; arc the densities of E[/ie„] and E[/i£;], respectively. 
In particular, when E is stationary (which implies that 0„ is also stationary), Ae„(2;) = 
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Ls„ e M+ and Xe{x) = Lb G R+ so that 

Note that, in general, 0„(a;) may be unbounded because the £'i(w)'s are not eompaet 
(e.g., Ei may be a random hne) or the union in not finite. Random elosed sets of this 
type are often taken to model real problems in applieations so that in the literature, 
many geometric processes like these are investigated, as point-, line-, segment- or plane 
processes, random mosaics, grain processes, etcetera (see, e.g., [3, 24]). In any case, it is 
commonly assumed in stochastic geometry that 0„ is locally finite (i.e., the number of Ei 
hitting any compact subset of R'' is finite a.s.); this is the reason why the compact window 
W CM.'^ is introduced in Theorem 4. Indeed, denoting by $^ the random number of E'i's 
hitting W and considering the restriction of 0„ to W, we may represent it as a finite 
union of random compact sets, almost surely, (and thus as in (20)) in the following way 



w 



where E^ = E,r\W. 

In the case of unbounded random sets, we consider the so-called Poisson line process^ 
as a simple example of applicability of Theorem 4. 

Example 3 (Poisson line process). Let Oi be the random closed set of Hausdorff 
dimension 1 a.s. associated with a Poisson line process in the plane (see, e.g., [24] for 
details). 0i is then a random collection of lines in such that the mean number of lines 
hitting any compact planar set is finite. Without any further requirement on Qi (e.g., 
stationarity or isotropy), it is easy to see that the hypotheses of Theorem 4 are satisfied 
with Y = V} (Oi n W^i). Indeed, by choosing 

HHei(c.)niyein-) ^ 
Hi(e,(.)ni^ei) ^''^ 

it follows that 

^^^Brix)) > ^^.^ , \^^^^ . r Vx e ei(c.) n W^ei,Vr e (0, 1), 
7^^(61(0;) n H/"ei) 

and, denoting by <I>^ the random number of lines hitting W^gi, 

E[H^(ei n W^ei)] < CEt**^] < 00, 

for a suitable constant C which depends on the dimension of the window W . 

Well known in the literature is the particular case in which Qi is stationary and 
isotropic. In this case, Qi can be described by a point process in the cylinder C* = 
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{(cosa,sina, p) -.p GM.,a G (0,7t]} in R^, with intensity measure A(d(p,a)) = L ■ dp ■ ^ 
for some constant L > 0. It is proved that the mean density Aei of 81 is equal to L and 
that the number Nr of lines of &i hitting the ball Br{0) is a Poisson random variable 
with mean 2rL (see [24], pages 249-250). We may note that the same result for the mean 
density can also be obtained as an application of Theorem 4: by (8), we have that 

^. P(SnB,(O)^0) ^. P(iV.>l) ,. l-e-2-^= ^ 

Ae, = hm = hm = Inn = L. 

rio 2r rio 2r rio 2r 

Let us observe that in the above example, in order to apply Theorem 4, for all w G i7, 
we have chosen the normalized measure rj as the restriction of the Hausdorff measure 
to 8i(w) n Wqi. Indeed, problems may arise in identifying a measure 77 needed for 
the application of the quoted theorem; in some cases, a proper choice of rj can be made 
by referring to another suitable random set which contains the relevant random set 6„. 
We further clarify this procedure by means of the following example, which suggests, in 
particular, a method to study more complex random sets, as briefly discussed at the end 
of the section. 



Example 4 (Segment process). Let 6i — [J^Si be a collection of i.i.d. random seg- 
ments Si in R'' with finite expected length, such that the mean number of segments 
hitting any compact set is finite. 

Given a compact window W C Mf^, if, for all a; G f2, we choose ry as in (28), then it is 
easy to check that Tw{&i) > mm{l,L}/n^{ei n W(si) =: l/Y, where 

L:= min {n^SA}. 

With respect to Theorem 4, the above is not a good choice for 77, because, unless L is 
bounded from below, we may well have E[Y] — 00; in this case, a possible solution to the 
problem is to extend all the segments having length less than 2 (the extension can be 
done homothctically from the center of the segment, so that measurability of the process 
is preserved). In particular, for any w € O, let 



Siiuj) 



S^iL0), iinHS^{uJ))>2, 

Si (uj) extended to length 2, if n^{Si (w)) < 2, 



and 81 := IJ^ Si. In this way. Theorem 4 applies with Y :— H^{Qi D W^i), by now choos- 
ing r] as in (28) and replacing 81 fl W^i with 81 fl W^i (see [25], Chapter 4 for details). 

Remark 11. A well-studied segment process in the literature is the so-called stationary 
Poisson segment process in (see [3, 24]). By proceeding as in the above example, the 
known result concerning the mean density of this kind of process can be rcobtaincd by 
applying (8) in Theorem 4. 

Example 4 provides a method of applicability of Theorem 4 to more complex random 
closed sets 8„ with Hausdorff dimension 1 < n < d a.s. (see also [25], Chapter 4 for an 
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application to Boolean models of spheres). An important class of 1-dimensional random 
sets are the so-called fiber processes (i.e., random collections of 1-rectifiable curves [3]). 
We recall that a set S CR"^ is said to be ri-rcctifiablc if it is representable as the image 
of a compact set K C K", with / : R" — >■ R'' Lipschitz, and we point out that condition 2 
is satisfied with ??(•)= ^"^{8 D •) for some closed set 5* 3 S' if / admits a Lipschitz inverse 
(see [1], page 111). Hence, if Oi is a sufficiently regular fiber process, an argument similar 
to that in the above example might be applied, by considering as 0i the random closed set 
given by the union of suitably extended fibers, so that information about the measure 
E[^ei] might also be obtained under hypotheses of inhomogeneity of the process. A 
relevant real system which can be modelled as a fiber process is the system of vessels 
in tumor-driven angiogencsis; estimation of the mean length intensity of such a system 
is useful for suggesting important methods of diagnosis and of dose response in clinical 
treatments [5, 8]. 
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